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We present a theoretical explanation for the interfacial zigzag instability that appears in 
anisotropic systems. Such an instability has been experimentally highlighted for an Ising wall formed 
in a nematic liquid crystal cell under homeotropic anchoring conditions. From an envelope equation, 
relevant close to the Freedericksz transition, we have derived an asymptotic equation describing the 
interface dynamics in the vicinity of its bifurcation. The asymptotic limit used accounts for a strong 
difference between two of the elastic constants. The model is characterized by a conservative order 
parameter which satisfies a Cahn-Hilliard equation. It provides a good qualitative understanding of 
the experiments. 



The zig-zag instability, undergone by straight rolls in 
two-dimensional extended systems like Rayleigh-Benard 
convection |l|] or electroconvection [|j fluid systems, have 
been extensively studied over the last decades. A similar 
zig-zag instability affecting anisotropic interfaces has re- 
cently aroused considerable interest in different fields, for 
example, in gas discharge system 0, crystal growth [||, 
rifts in spreading wax layer [|| , Ising wall in nematic and 
cholesteric liquid crystals and chevrons layer structure 
of smectic liquid crystals |j. In this letter, we present 
a theoretical explanation for this interfacial zigzag insta- 
bility and show how it justifies the Ising wall dynamics in 
liquid crystal samples. Besides, the instability is experi- 
mentally characterized for a thin layer of nematic liquid 
crystal homeotropically anchored. From an amplitude 
equation, which is valid in the vicinity of the Freedericksz 
transition, we derive a solvability equation that describes 
the interface dynamics at the onset of the bifurcation. 
This model is in qualitative agreement with the experi- 
mental results, even though these latter were obtained far 
from the Freedericksz transition. Numerical simulations 
allow to strengthen the relevance of the model. 

In two-dimensional extended systems, the dynamics 
of a straight interface can be entirely characterized by 
the dynamics of the interface position P, function that 
parametrizes the interface in space and time. Henceforth 
we consider that the variables x and y describe respec- 
tively the transversal and longitudinal directions of the 
straight interface. As a consequence of the translational 
invariance of the interface (P — ^ P -h Po) and of the 
space reflection symmetry in the tangential direction to 
the interface (\j — > — y), the order parameter P{t^y) 
shows a spatio-temporal evolution which is modelled by 
the Kuramoto-Sivashinky equation [||, as long as one as- 
sumes that, the interface only depends on the tangential 
direction. If one considers moreover the space reflection 
symmetry in the transversal direction (P — > — P ), then 
the order parameter satisfies the following diffusion equa- 
tion 



When s is small (positive or negative), the order param- 
eter dynamics is described by the asymptotic equation 



dtP = ePy 



Pyy 



(2) 



where £ is the diffusion (antidiffusion, £ < 0), PyPyy the 
nonlinear diffusion and the last term the hyper diffusion. 

In the preceding expressions, we have only considered 
transversal perturbations to the straight interface and 
adiabatically eliminated the diffusive behavior in the lon- 
gitudinal direction, approximations that are both correct 
for anisotropic systems. In isotropic systems, a straight 
interface can be linearly unstable {e < 1) along the lon- 
gitudinal direction as a result of the interface symmetry. 
However, in this case, nonlinear terms appear simulta- 
neously {dxPdyyP, and dyPdxyP evaluated at the in- 
terface) that destabilize the solutions of eq. (||). Due 
to these nonlinear terms, the vertices, shown on figures 
^ and ^, will become fingers that propagate along the 
transversal direction [^. Therefore all physical systems 
that show a zig-zag instability leading to a zigzag coars- 
ening dynamics are anisotropic [^, H, 0] • 

In eq. (^) we have assumed that the nonlinear and 
fourth-order derivative terms saturate the instability, 
which imposed the signs of these two terms. Note that 
this equation is a continuity equation that expresses the 
conservation of the area enclosed by the curve P{y^ t) and 
the y-axis. This law stands from the fact that the inter- 
face connects two symmetrical states (P — > — P) and 
hence, when one part of the interface moves in one direc- 
tion, another part moves simultaneously in the opposite 
direction. Moreover, this equation is variational, i.e. it 
can be rewritten in the following form 
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(3) 



where the free energy T only depends on the order pa- 
rameter derivatives. Introducing the variable A 



y^ 



dtP = ePy 



(1) 



which describes the local tilt of the interface with re- 
spect to the original orientation of the straight interface, 
the latter equation is reduced to 
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This equation is the weh-known Cahn-Hihiard equation 
[|l2|, which describes the phase separation dynamics in 
conservative systems. The dynamical behavior of this 
equation has been studied by many authors (see for ex- 
ample [|l^, Q and the references therein). It is notewor- 
thy that the relevant conservation law for the dynamics 
is actually 



M = J A{y,t) dy 



(5) 



and not the integral of P, as a consequence of the trans- 
lational invar iance. When M is zero, the dynamical evo- 
lution, deduced from the theoretical analysis, is in good 
quantitative agreement with the experiments, as we shall 
see later. A wave number of order ^/—e/2 arises initially 
in the system. Later on, the system exhibits a coarsen- 
ing dynamics (see fig. ||) during which the wavelength 
increases. New periodic solutions emerge with a loga- 
rithmically slow growth rate 

To find the steady solutions, one has to minimize the 
free energy taking into account the area conservation. 
This can be done by considering a Lagrange multiplier 
(A). Hence, the free energy takes the form 



dy ■ 



A2 



A* 
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(6) 



The family of solutions that minimize this free energy 
satisfies the differential equation 



A^^ = A + eA + A^ 



(7) 



Obviously such solutions also satisfy eq. (H). The global 
minimum of the free energy Q is either the homoclinic 
solution (which is an interface with three facets) or the 
constant solution. This one, as well as the facets, can 
be any of the local minima of the potential V = AA + 
e^ + ^ . We remark that these solutions have different 
Lagrange multiplier values. 

Let's consider now an anisotropic nematic liquid crys- 
tal sample that exhibits experimentally an interfacial zig- 
zag instability, as emphasized below. This instability will 
be theoretically described in the framework of the non- 
linear elasticity theory of liquid crystals. 

In the experiments, we used a cyanobiphenyl com- 
pound (5CB) whose anisotropic physical properties at 
25°C are: elastic constants Ki = 6.3, = 4.1, = 8-4 
(10~^ dynes); dielectric anisotropy = 11.3; diamag- 
netic anisotropy Xa — 1-142 (10~^cm^^~^); rotational 
viscosity 71 ~ 10~^Pa.s. 





FIG. 1: Experimental set-up: the sample of liquid crystal is 
placed between two permanent magnets which determine the 
homogeneous magnetic field H. The physical phenomena are 
observed through a polarizing microscope. 



The samples are made of two glass plates separated 
by thin mylar spacers that determine the cell thick- 
ness (between 50 and 250 /im). The glass surfaces are 
treated with lecithin to provide homeotropic anchoring. 
The sample is subjected to a sinusoidal vertical electric 
field E — Ecz {Veff. ~ — 9^) with a high frequency 
(^ bkHz) in order to avoid charge injection or electrocon- 
vection phenomena. Moreover, two permanent magnets 
induce a homogeneous horizontal magnetic field in the 
sample. The field magnitude can be changed by bring- 
ing the magnets nearer or farther. Its maximum value is 
about 0.55 Tesla. 

The application of a horizontal magnetic field can in- 
duce, for a high enough magnitude, a partial reorien- 
tation of the bulk molecules along the field direction. 
Due to the twofold degeneracy of the bifurcated state, 
domains of opposite orientation may be created [0. 
The interface between a pair of domains is called Ising 
wall [|ll|. An Ising wall in a splay-bend configuration 
is formed by using the fiux lines curvature in a region 
where the field is inhomogeneous (see fig. |l|). Then, it is 
quenched into the area between the two magnets, where 
the field is homogeneous. Its width is determined by the 
distance to the threshold of the Freedericksz transition 
and can be modified by variations of the fields magni- 
tude. 

The interface behavior is observed through a polarizing 
microscope (see fig. |l|). Video films or numerical images 
(see fig. ^ can be registered thanks to a 3CCD camera 
placed on the top of the microscope. 

The experiments have been carried out far from the 
Freedericksz threshold, since near the transition any little 
imperfection on the parallelism between the sample and 
the magnetic field plane makes the wall drift towards one 
edge of the cell. 
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FIG. 2: Spinodal decomposition of the interface observed 
through a microscope with crossed polarizers. 



The splay-bend Ising wall is ever unstable, for the con- 
trol parameters that we used experimentally, when it is 
thrust into the homogeneous area of the magnetic field. 
Initially, the interface develops an instability character- 
ized by a well-defined wavelength, which is determined by 
the experimental parameters. Later on, the sinusoidal in- 
terface becomes an angled line composed of pieces of wall 
turned with angles (see fig. ^. Two adjacent pieces, 
whose orientations are opposite, are connected by a re- 
gion of strong curvature of the line that we called kink. 
The dynamics consists then in reassembling domains of 
even orientation, the angle of the "zig" and "zag" facets 
staying unchanged (coarsening dynamics). This process 
occurs thanks to annihilations of kinks and without char- 
acteristic lengthscale. Actually, the averaged domain size 
increases regularly in time (see fig. The dynamics, 
which tends to separate the zig and zag facets, is the 
one-dimensional counterpart of the spinodal decomposi- 
tion dynamics observed in conservative binary mixtures 

The instability is triggered by the elastic anisotropy of 
the liquid crystal whose influence is emphasized in the 
wall. Indeed, the involved distortions depend on the ori- 
entation of the interface with respect to the magnetic 
field direction. For most of the usual compounds, the 
energy cost of a wall aligned with the field (twist wall) 
is lower than the one of a wall perpendicular to the field 
(splay-bend wall) since K2 < Ki . Consequently, the sys- 
tem (splay-bend wall) reduces its energy by changing the 
direction of the wall. A global rotation of the interface 
can not occur in an infinite medium. Therefore, the inter- 
face is forced to rotate locally and is divided into facets 
turned with opposite angles ±^o- These angles, whose 
values are determined by the control parameters (see fig. 
H), result from several effects. More precisely, the rota- 
tion is favored by the elasticity but is made difficult by 
the resulting elongation of the interface along its orig- 
inal direction and by the escape of the molecules from 
the vertical plane containing the magnetic field. The lo- 
cal reorientation generates many defects (kinks), which 
will eventually disappear since this leads to a decrease in 
energy. 

The experiments have been performed far from the 



FIG. 3: Experimental curve giving the slope Py of the in- 
terface versus the dimensionless ratio 7^ is the electro- 
magnetic coherence length and the magnetic coherence 
length) . 



Freedericksz transition. However, for the sake of simplic- 
ity, the theoretical approach has been developed near this 
transition. At the onset of the bifurcation, the system is 
described by the two-dimensional Landau equation 



drZ = eZ- bZ^ + {Kidl 



- K^d^) Z. 



(8) 



where £ = -CaE^ -Ks^^XaH^ . ^ = 

4 



1 ~ 2^3; 

\eaE'^ + ^XaH^, T = t/71, and the director takes the 



form n = (^Z cos(^), 0, 1 - ^ cos2(^)y This equa- 
tion contains wall-type solutions. When both diffusion 
constants {Ki/ji and i^2/7i) are of the same order of 
magnitude, these solutions are always stable. In order 
to explain the appearance of the interfacial instability, 
one has to consider that K2 is much smaller than the 
other constants Ki and Ks {K2 ~ e). Then, the solution 
Z = is marginal with respect to perturbations in the 
^/-direction. Note that the interface identifies with the 
center of the Ising wall, which is a curve where Z = 
(homeotropic state). Hence, higher-order terms must be 
considered in order to study the stability of the interface, 
and the dynamical equation reads now 



drZ = sZ- bZ^ + {Kid'^ + K2dl) Z - 



-^d,,y.Z^^Ks [Z{dyZf-^dyyZp) 



a 

Z^ 
2 



d^lylZ 



where we have considered the asymptotic limit Z ~ e^/^, 
dlZ - dlZ - £3/2, Ki Ks ^ 1, and K2 - e, 

with £ < /i < 1, a =eaE^ + Ks^ and ^ = (i^ - 2^). 
Then, the expression of the local director is given by 




V 



■ cos^ 



In this asymptotic limit, the Ising wall solution reads 
Z = Y^^/^ tanh (^■\/e/2Ki {x — Xq)^ • In order to describe 
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FIG. 4: Numerical simulation of the eq.(8) for a high 
anisotropy of elasticity : ^ = 0.056 , K3 = 1 , a = 1. The 
picture shows the Y-component of the director in the system. 



the role of the elastic anisotropy. Such a limit contains all 
the ingredients required to give an adequate description 
of the experimental observations. 
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the dynamical behaviour of the wall under perturbations, 
we introduce the following ansatz 



Z {x, P{y, T)) = ^^ tanh (^^ - P{y. T)j + w{x, P) 

(10) 

where w{x^P) is a small perturbative term. One can 
show that, P satisfies the following solvability equation 



OtP = DiPy 



D2PyPyy 



D.P 



2,-Lyyyy 



(11) 
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2Kis 
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with Di ; 

Thus, the characteristic time-scale of the interface dy- 
namics (/i~^6:~^) is smaller than the time-scale of the 
interface formation {s~^)^ which is consistent with ex- 
perimental observations. The expression of Di indicates 
that the elastic anisotropy can give rise to the instability 
(1^1 < 0). Numerical simulations of eq. (^) confirms the 
latter calculations (see fig. ^). 

As a consequence of the finite size of liquid crystal 
sample, one experimentally observes a global rotation of 
the interface, which is superimposed to the coarsening 
dynamics. This favours one orientation of the domains, 
but does not change the angle of the facets. 

In conclusion, we have presented a theoretical expla- 
nation of the interfacial zig-zag instability that appears 
in anisotropic systems. The instability of an Ising wall 
formed in a nematic liquid crystal has been more specifi- 
cally investigated, both experimentally and theoretically. 
Close to the Freedericksz transition, this instability is de- 
scribed by the Cahn-Hilliard equation, and leads to the 
facets separation and coarsening dynamics. The asymp- 
totic limit used in the derivation of this model emphasizes 
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